Confronting quintessence models with recent high-redshift supernovae data 
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We confront the predictions of different quintessence models witfi recent measurements of the 
luminosity distance from two sets of supernovae type la. In particular, we consider the 157 SNe la 
in the Gold dataset with z < 1.7, and the more recent data containing 71 supernovae obtained by 
the Supernova Legacy Survey (SNLS) with z < 1. We numerically solve the evolution equations 
for the Ratra-Peebles inverse power law model, the double exponential and the hyperbolic cosine 
quintessence models. We obtain confidence regions from the two datasets in the Hm—ci and Qm—w^ 
planes for the different models and compare their predictions with dark energy models with constant 
equation of state. 

PACS numbers: 98.80.-k, 98.80.Es 



INTRODUCTION 



The growing observational evidence from high-redshift 
supernovae 0, Q and other cosmological data Q] sug- 
gests that the dominant component of the universe today 
is some sort of dark energy fluid with negative pressure 
[J, |5j- Indeed a component with constant equation of 
state px — wxpx and wx < —0.78 [3j fits the exist- 
ing data with reasonable goodness (XCDM model). In 
particular, it includes the most economical explanation 
for the present state of accelerated expansion of the uni- 
verse, i.e., the presence of a pure cosmological constant 
with Wx = — 1- However, ther are other models in which 
the equation of state depends on redshift and which are 
also able to fit the data with comparable quality. This 
is the case of the so called quintessence models [a, H , in 
which it is the presence of an evolving scalar field with 
appropriate potential term what plays the role of dark 
energy. 

Quintessence models exhibit an interesting property 
usually called tracking behaviour |3|. This means that 
there is a wide range of initial conditions for which the 
evolution of the scalar fields converges to a common evo- 
lutionary track. Furthermore in such a tracking regime 
the equation of state of quintessence w^ remains almost 
constant. Moreover there are particular models, usually 
called scaling models |9|, llOl lll| , in which w^ mimics the 
equation of state of the dominant component through- 
out the cosmological evolution. These properties allow 
quintessence models to alleviate the fine tuning problem 
of XCDM models, although in any case the scale of the 
quintessence potential has always to be tuned in order to 
reproduce the data. 

Apart from quintessence, other models have been pro- 
posed in the literature. Thus for instance, scalar fields 
with non-canonical kinetic terms (k-essence) 12]; the 
generalized Chaplygin gas model [l3|, which in principle 



allows for an unification of dark energy and dark mat- 
ter; infrared modifications of General Relativity, as for 
instance in extra dimensional theories [1J| or the modi- 
fication of the Friedmann equation in the so called Car- 
dassian models jl^, lla] ■ 

The present SNe la data are not very constraining 
when trying to determine the nature of dark energy and 
different dataset favor different regions of the parame- 
ter space. Thus for instance, it is known that the Gold 
dataset Q prefers models with wx < — 1 [l3>ll3l) the so 
called phantom dark energy models, whereas for the more 
recent Supernova Legacy Survey (SNLS) l3 , the best fit 
is closer to a pure cosmological constant []j|, |2flj . There- 
fore it is worth exploring what kind of constraints are 
imposed in each case, and this is the aim of the present 
work. However the information we will obtain will be 
very limited because even the highest redshift points 
in those sets explore only relatively recent epochs, and 
therefore the present supernovae data alone are not able 
to discriminate between the different models (although 
larger future supernovae catalogues are expected to im- 
prove in an important way the present constraints). Ac- 
cordingly we should rely on complementary information 
in order to distinguish quintessence models from a cosmo- 
logical constant or other alternatives |2lU22l |. Thus for in- 
stance, in the particular case of quintessence, dark energy 
is generated by a scalar field, and there exists the possi- 
bility of having density or velocity perturbations which 
could affect CMB anisotropics [231 . In constrast, the pure 
cosmological constant case does not support such pertur- 
bations. 

In this work we will concentrate on different 
quintessence models and derive constraints on their 
parameters from the two mentioned Gold and SNLS 
datasets. We consider three models: an inverse power- 
law Ratra-Peebles (RP) model 0, V{cj)) = M^+°'/cj)°', the 
double exponential (DE) potential y((/)) = M'*(e"'^-Fe^"^) 



|2 4L |25j which exhibits scaling behaviour and the hyper- 
boHc cosine (HC) V{(j)) = M4(cosh(A(?!)) - 1)" jMJ, which 
posseses an oscillatory behaviour at recent epochs. Un- 
like previous works |27l |28J , we do not parametrize the 
equation of state ^^(z) in order to obtain explicit expres- 
sions for the luminosity distance-redshift relation dhiz). 
Instead, we solve numerically the evolution equation for 
the scalar field and the universe scale factor. This allows 
us to numerically obtain d^ iz) for each model. The value 
of the fitted cosmological parameters is known to have 
a strong dependence on the particular parametrization 
chosen for the equation of state (see [3, |29| and refer- 
ences therein). In particular, at least three parameters 
are needed in order to take into account properly all the 
information encoded in the SNIa data. Our approach 
does not rely on particular parametrizations, but instead 
we consider the full redshift dependence by numerically 
solving the evolution equations. 

The plan of the paper goes as follows, in section 2 we 
review the models properties and obtain the appropri- 
ate equations of motion. In section 3, we compute the 
luminosity distance expressions and perform the corre- 
sponding cosmological fits. Section 4 contains the main 
results of the paper with the confidence regions for dif- 
ferent spaces of parameters and finally we include a brief 
section with conclusions. 



QUINTESSENCE MODELS AND EVOLUTION 
EQUATIONS 

Inverse power law potential 

In the first model that we will consider, the potential 
reads: 



V{^) 



^^4+Q 



(1) 



This model has a tracker attractor, but it is not scaling, 
although the equation of state is determined by that of 
the dominant component: 
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where, for a > —2, in the radiation or matter eras we 
have w^ < wr^m respectively, i.e. the energy density de- 
creases more slowly than the dominant component and 
it eventually becomes dominant at late times. The cor- 
responding equations of motion reduce to the Friedmann 
equation for the scale factor, together with the scalar field 
evolution equation: 



H' = 




V{(j))+ J2 P^^oc 



-3{w, + l) 



i=M,R 



4> + 3H^ + V'{(I)) =0 



(3) 



where pi^/ po = i~li with i = M, R and po is the critical 
density. Here a dot denotes derivative with respect to 
cosmological time. 

In order to solve these equations numerically, we trans- 
form them into dimensionless equations, defining: r = 
Hf)t and (j) — (87rG)^/^(/). Now a prime will denote deriva- 
tive with respect to t: 
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+ V{4>) + VtRa-^ + VImc 



(4) 



where 



y^^) = 4 



A = 






a+2 



(5) 



The initial conditions for (j) and a will be given by: (t>(0) = 
5 • 10~\ 4>'{0) = and a(0) = IQ-^. We wiU explore the 
parameters range a — — 5.5 and A = 0.3 — 8-10^. Notice 
that for those values of the parameters, the initial scalar 
field energy density is a small fraction of pm and pR. 
Notice also that thanks to the attractor behaviour, the 
late time evolution is not very sensitive to the particular 
initial conditions chosen. 



Double exponential potential 



In this case: 



V{4>) = M*(e"'"*" + e'""^) 



(6) 



where k = vovrG. The energy density in this model 
follows a scaling behaviour at early time, whereas it be- 
comes dominant at late times with appropriate equation 
of state: w^ ~ —1. The corresponding dimensionless 
potential now reads: 



F(0) 
A 
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(7) 



and the corresponding equations of motion are given in 
(gl). The initial conditions are 0(0) = 0.5, 0'(O) = 
and a(0) = 10^"^. The range of parameters considered 
are /3 = 20, a = 0.1 - 1.5 and A = QA- 120, where 
for simplicity we have fixed the constant (3. The a and 
(3 constants essentially fix the value of w^ today, and 
we have checked that fixing (3 still allows us to cover a 
wide range of values in w^, just varying a. Again for 
those particular values, the initial energy density in the 
scalar field is a small fraction of the radiation and matter 
densities. 



Hyperbolic cosine potential 



In this case the dimensionless potential reads: 
F(^) = A{cosh{X4>) ~ 1)" 



ml 



(8) 



In the limit A^ ^ 1, i.e. at early times, the potential be- 
haves as a single exponential and it is possible to choose 
the parameters so that the model posseses a scaling be- 
haviour. In the opposite limit, \4> <C 1 (late times), the 
potential can be approximated by Vi^cj)) ex (A0)^", and 
the scalar field oscillates around (/) = 0. When oscilla- 
tions start, the equation of state no longer mimics the 
dominant component, but it also start oscillating with 
average value given by: 



(W0> 



a-l 
a + 1 



(9) 



Accordingly, the value of a determines the present equa- 
tion of state. For a = 1 the oscillations behaves as nonrel- 
ativistic matter, whereas if we require accelerated expan- 
sion today, i.e. w^ < —1/3 then we should have a < 1/2. 

The numerical integration of equations Q now has an 
additional difficulty. In the interesting case (a < 1/2), 
the potential derivative V appearing in the first equation 
diverges at the origin. Since both and (/)' are continuous 
at the origin, in order to avoid numerical instabilities we 
have smoothed the divergence modifiying the potential 
as y(0) = A(cosh(A0) - (1 - e))". We have checked that 
for sufficiently small values, the results do not depend on 
the e parameter. 

The initial conditions considered are (j){0) — 1.6, 
^'(0) = 0, a(0) = 10~^. The free parameters are 
a = 0.02 - 0.4, A = 0.3 - 14 and we have fixed Xa = 5, 
which ensures again that initially quintessence is not the 
dominant component of the energy density and that it 
follows a scaling behaviour. 



COSMOLOGICAL DATA FIT 

After numerically solving Q , we obtain a discrete time 
evolution for the scale factor a{t) and the scalar field (j){t) 
for a given set of parameters {A, a). The time depen- 
dence will be used to compute the theoretical values of 
the luminosity distance to each SNe la with a given red- 
shift z. Thus we use the well known expression, derived 
from the FRW metric under flat prior, expressed in terms 
of the dimensionless time variable t. 



dL{T) 



dr 



a(n) Jro a('^) 



(10) 



The observations of supernovae measure essentially the 
distance modulus /Lt, which is the diference between the 



apparent magnitude m and the absolute magnitude M, 
and relates to the luminosity distance as: 

^ith = m-M = 51ogdL + 51og(^^^^W25 

= 5logdL + M (11) 

The M value can be assumed constant once the necessary 
corrections are applied on m{z). 

Assuming tq to be the value for which a(ro) = 1, the 
scale factor depends on redshift as a{T) = (z + 1)^^. 
As a consequence, for each supernovae in a given set, 
we derive a(Ti) from its redshift and obtain ti from the 
numerical output data of Ql . Hence the integral H10|) can 
be numerically evaluated to obtain fith for each {A, a), 
once the value of Hq is fixed. 

Once we obtain fith, the comparison to its observa- 
tional value will enable us to carry out a x^ statistical 
analysis. For this purpose, we have considered two sets 
of supernovae. On one hand, the Gold set, compiled by 
Riess et. al. [l|, containing 143 points from previously 
published data, plus 14 points with z > 1 discovered 
with the HST, all reduced under the same criteria in or- 
der to improve the errors arising from systematics. On 
the other, the SNLS set, comprising 71 distant super- 
novae (0.15 < z < 1) discovered during the first year of 
the Supernova Legacy Project (SNLS) 2], alongside with 
44 SNe la from other sources that feeds the nearby zone 
(0.0015 < z < 0.125), and which are also included in the 
Gold set. 

The process followed to obtain x^ is slightly different 
for each set. In the Gold set we have used the obser- 
vational distance modulus fj-exp given by Riess et. al. 
together with its associated error cr^ to compute x^ as a 
function of the free parameters of the quintessence model. 



x'(Aa,M) 
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(Mohs - ^J.thiA, a, M))"^ 



(12) 



The dependence on Hq has been accounted for in l(T^ 
through the nuissance parameter M , which is indepen- 
dent of the data points and the data set. We have 
marginalized x^ over all values of M by ex panding and 
minimizing ifT^ with respect to M (see [la. Il7lll9|| 1 . 

On the other hand, for the SNLS data, a more detailed 
relation between ^ and the observational measurements 
has to be considered in the statistical analysis, which 
implies recursively fitting two new non-cosmological pa- 
rameters in the calculation of x^. 

The expression for the observational distance modulus 
used by the SNLS team includes these two parameters to 
measure the impact of the rest frame color parameter (c) , 
and the light curve strecht (s), on the distance modulus. 



I^'oil'' = -(-) 



M + ai{s - I) - a2C (13) 



Introducing this equation into the expression for x^ 

x'(a,^,ai,a2,Af + M)= (14) 



EN 



we obtain a six parameter dependence (five parameters in 
practice, since we can replace M + Af by a single additive 
constant) We minimize following the process suggested in 
[21 , i.e. marginalizing with respect to ai , 02 and M + M 
for each pair of (A, a) values. Notice that marginalizing 
over M or M + M is equivalent to marginalize with re- 
spect to Hq. Again /i and cr^ are obtained from the data 
in0. 



RESULTS AND MODEL COMPARISON 

The best resulting cosmologies obtained after fitting 
each model to the Gold and SNLS sets are summarized 
in Table I and II respectively. In each sample, we find 
the same x^ for the best fit with the three quintessence 
models considered in the paper, i.e. Xmin = 177.0 for 
the 157 SNe la of the Gold set, and xLin = Hl-O for 
the 115 SNe la of the SNLS set. In addition, the min- 
imun x^ foi' all models corresponds to a = 0. For this 
value the RP and HC potentials turn into a cosmological 
constant term and we have used this fact to check our 
results, finding good agreement for the Qm value when 
compared to the Riess et. al., and Astier et. al. XCDM 
model. On the other hand, albeit the DE potential does 
not strictly behave as a A term, due to the nonvanish- 
ing exponential, the uj^ value tends asymptotically to -1, 
leading to the same result as DE and HC. This implies 
that a quintessence potential, and a pure ACDM model, 
can not be distinguished only by fitting {Q]\i, i^ip), since a 
quintessence model can be tuned to accurately resemble 
a cosmological constant, at least in the redshift interval 
explored by the actual SNe la sets. In any case, notice 
that prj, evolves in time for any non-zero value of a for 
the three potentials. 





Gold 


1 




^M 


i^4> 


a 


x' 


RP 

HC 
DE 


0.30±^-i'? 

o.30±g:;g 
o.30±EJ:8? 


-l±"l^^ 
_^^0^4 

_^^0^22 


0^1.35 

0^1^08 


177.0 
177.0 
177.0 



TABLE I: Best fitting cosmological parameters, with their 
associated 68% errors for the Gold set. The degeneracy in a = 
leads to the same value of Q.m for the three potentials in a 
set. Moreover, due to the form of the quintessence Lagrangian 
values of uj^ below -1 cannot be obtained. 

Combining the data obtained by solving @ for each 
potential, we may relate (A, a), and its corresponding 
fitting x^; to a single (w^jJIm) pair which allows us to 



plot the confidence regions for the two parameter combi- 
nations {a,Q.M) and {uJcjj,Q.m)- 

For the HC model, we use directly (oj^) = j-?, in- 
stead of the value derived from eq.Q, due to the dif- 
ficulties arising from the calculation of (w^)""™ in the 
cases when the oscillations fails to complete a whole pe- 
riod. Moreover, a comparison between {oj^) and (w^)""™ 
shows a small average discrepancy, Aw^ = 10"'^, which 
justifies our approximation, and serves as a check of the 
numerical solution. 



RP 
HC 
DE 



SNLS 



Q. 



M 



<0.36 
<0.36 
<0.36 



i^,)> 



<-0.58 
<-0.61 
<-0.54 



>0 
<0.3 
<1.27 



X 



111.0 
111.0 
111.0 



TABLE II: Best fitting cosmological parameters, and upper 
bounds at 95% C.L. for the SNLS set. The constraints on the 
cosmological paramaters hold for Q,m > 0.06 for the RP and 
DE potentials, and ^m > 0.12 for the HC. 

Fig.l to Fig. 3 show the confidence regions for the mod- 
els in the a — Q.m and w^ — Qm planes, together with the 
XCDM contour plots. The XCDM regions have been cal- 
culated using the same code on the Friedmann equation 
with a dark energy term {px,^x) instead of a coupled 
scalar field. In the a— Om plots the contours grow asymp- 
totically for decreasing VLm JS^]- Nevertheless, as uj^ fol- 
lows the same growing trend, it is possible to restrict the 
maximum value of ujip . 

At the 95% confidence level and for VLm > 0.1, we find 
,9oid / n r;R ,,,SNLS / _o_58^ for the RP model, and 

-0.56, for the DE model. The 



qold 



< -0.56, c^f^^^ 

< -0.55 ''^^^^ 



< 



i^A 



< 



HC potential also shows similar results for J^m > 0.15, 
limiting w^°''' < -0.58, and lu^'^^''^ < -0.64. Notice that 
we do not extract any result for values below JIm =0.1 for 
RP and DE or fluj =0.15 for HC, since the computational 
effort increases in those regions. In fact in the HC model 
that region is also highly sensitive to the value of the 
smoothing factor, e. 

By comparing the plots for the two datasets, we 
see that the width of the contours is essentially the 
same. However, the SNLS 68% C.L. contour never closes 
for Um >0.06, whereas the Gold set constraints the 
cosmological parameters, for all models, within 0.25< 
i^M <0.36 and -1< uj^ <-0.75, at this significance level. 
For the RP potential, these results can be compared to 
analysis done with previous supernovae data in |2ll l80l | . 

Concerning the XCDM contours, one remarkable fea- 
ture is that as w^ 1-^ —1, the SNLS confidence regions 
for quintessence and XCDM tends to overlap. This is an 
expected fact, since the best fitting value for XCDM ac- 
cording to Astier et. al., is uj^ =-1.02, quite close to the 
quintessence best fit, which is precisely the degenerate 
case a = 0. On the other hand, the best fitting value 
for the Gold set, assuming no priors, is w^ =-2.3. This 



implies an important difference in ^ with respect to the 
quintessence models, which prevents the confidence re- 
gions from matching. In spite of this, it is noticeable the 
similarity between contours given the fact that the 68% 
C.L. region of XCDM is out of the plotted interval. Fi- 
nally, for all models, the quintessence 99% contours place 
an outer boundary containing all the XCDM regions, im- 
plying that any of the XCDM cosmologies, might be ob- 
tained from a quintessence potential. 



CONCLUSIONS 

We have compared the predictions of three 
quintessence models with the measurements of the 
luminosity distance vs. redshift from two SNe la 
catalogues, (Gold and SNLS surveys). We have obtained 
the corresponding confidence regions in the {ol^^m) and 
(w^jJIm) planes and compared them with the predic- 
tions of dark energy models with constant equation of 
state (XCDM). From the plots we see that the different 
confidence regions for quintessence and XCDM models 
cover essentially the same areas in the (flM,^^) plots. 
This makes it difficult to discriminate between the 
different models just from SNe la data alone. In all the 
three cases, the best fit cosmology corresponds to the 
a = case, in which the quintessence potential reduces 
effectively to a cosmological constant. Although SNLS 
data favor a pure cosmological constant, they yield larger 
confidence regions than the Gold dataset. Bounds on 
the cosmological parameters iljv/, ^4> and the potential 
paramaters a have also been obtained. In principle it 
could be interesting to use the differences in constraints 
between the two sets as an estimate of possible residual 
systematics in the SNe la data. However, in practice, 
since both datasets have been obtained and reduced 
in different ways and possibly the actual data is still 
dominated by statistical errors, it would be difficult to 
perform such estimations. In any case this is beyond the 
scope of the present work. 
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FIG. 1: Confidence regions at the 68%, 95% and 99% CL. (corresponding to Ax^=2.3, 6 and 11.8 for a two-parameter fit) 
for {a,nM) (left panels) and ((jj^,f2A/) (right panels), in the Ratra-Peebles quintessence potential, obtained using the Gold 
(top panels) and the SNLS (bottom panels) samples. The dotted lines, on the right panels, show the probability contours for 
a XCDM model: at the 99% and 95% for the Gold set, and at the 99%, 95%, and 68% for the SNLS set. 
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FIG. 2; Confidence regions at the 68%, 95% and 99% C.L. (corresponding to Ax =2.3, 6 and 11.8 for a two-parameter fit) for 
(Q:,r2jv/) (left panels) and (aj^,fiAf) (right panels), in the Double Exponential quintessence potential, obtained using the Gold 
(top panels) and the SNLS (bottom panels) samples. The dotted lines, on the right panels, show the probability contours for 
a XCDM model: at the 99% and 95% for the Gold set, and at the 99%, 95%, and 68% for the SNLS set. 
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FIG. 3: Confidence regions at the 68%, 95% and 99% C.L. (corresponding to Ax^~2.3, 6 and 11.8 for a two-parameter fit) 
for {a,nM) (left panels) {uj^,flM) (right panels), in the Hyperbolic Cosine quintessence potential, obtained using the Gold 
(top panels) and the SNLS (bottom panels) samples. The dotted lines, on the right panels, show the probability contours for 
a XCDM model: at the 99% and 95% for the Gold set, and at the 99%, 95%, and 68% for the SNLS set. 



